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Abstract 

We derive the minimum observing time scales to detect a blazar at a given flux level 
with the LAT on GLAST in the scanning and pointing modes. Based upon Phase 
1 observations with EGRET, we predict the GLAST detection rate of blazar flares 
at different flux levels. With some uncertainty given the poor statistics of bright 
blazars, we predict that a blazar flare with integral flux £ 200 x 10~ 8 ph(> 100 
MeV) cm" 2 s _1 , which are the best candidates for Target of Opportunity pointings 
and extensive temporal and spectral studies, should occur every few days. 
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1 Introduction 



The statistical significance for detection of a known source at the na level is 
n = S/ vB, where S is the number of source counts, B is the number of back- 
round counts, and the diffuse background is assumed to be precisely known 
lj). The observer is also assumed to follow a predetermined data analysis 
protocol and not to make multiple rebinnings of the data. 

In this paper, we derive the minimum observing time for the Large Area 
Telescope (LAT) on GLAST (Gamma ray Large Area Space Telescope) to 
detect high-latitude 7-ray sources at given flux levels, and estimate the rate 
at which blazar flares at certain flux levels are expected to be observed with 
GLAST. 



In Sections 2 and 3, we present calculations of signal and background counts, 
respectively. The five count and the five sigma requirements for detection are 
derived in Section 4. EGRET blazar studies, described in Section 5, are used 
to infer the probable rate of flaring emissions that GLAST will measure, as 
discussed in Section 6. 
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2 Source Counts S 



The point-spread function (psf) of the LAT on GLAST is characterized by 
angle 6 psftG = (3.5°/57.3°) (E/E 100 )~ 2 / 3 = 0.061u- 2 / 3 , where E 100 = 100 
MeV and £7 = E(MeV) 0). We assume an energy- dependent analysis of the 
GLAST data, where the number of counts is integrated in the direction of 
a catalogued source over the solid angle AQ G = tt0 2 s/jG = 0.012m" 4 / 3 . For 
this analysis procedure, the fraction / 7 of detected source photons is therefore 
independent of energy. For an azimuthally-symmetric, Gaussian distribution 
oc exp[— (0/0 ps f,G) 2 } of count directions, / 7 = 1 — exp(— 1) = 0.63, though 
fry = 0.68 would apply if 9 ps f is defined as the 68% containment angle. 

The number of source photon counts with energy > E within AQ G of a given 
point source is, in this approximation,given by 

At oo 

S = f,JdtJ dE' A[E', 6(t), 0(t)] 4> S {E', t) (1) 

E 

where <p s (E,t) is the source photon flux (ph cm -2 s _1 E~ 4 ). We approximate 
the energy- and angle-dependent effective area of the LAT by the function 

A(E) = A u(6)(E/E 100 ) aie) = XA Q u ao , (2) 



where u(6 = 0°) = 1, a{9 = 0°) = do, and 6 is the angle between the direction 
of a source and the axis of the LAT, and <fi gives the azimuthal angle of 
the source with respect to the LAT orientation. Here we assume azimuthal 
symmetry of the detector effective area. An approximation for the effective 
area that satisfies the GLAST Science Requirements Document is Aq pa 6200 
cm 2 and a ~ 0.16 for 100 MeV ^ E £ 10 GeV (compare Foldout B in the 
GLAST proposal; this approximation breaks down above 10 GeV). For 20 MeV 
£ E ^ 100 MeV, a ~ 0.37. The net observing time is denoted by At = 10 4 t 4 
s, and the exposure factor X accounts for the fraction of time that a given 
source is being detected by the LAT. The effective area of the LAT drops 
by a factor-of-two at 9 — 55°, so that we take X = 0.2A .2 with X .2 — 1 
in the scanning mode, noting that |(1 — cos 55°) =21%. (Signal detection on 
a 10 4 second time scale is of special interest in blazar studies, because the 
light-travel time across the Schwarzschild radius of a 1O 9 M is ~ 10 4 s.) 

Normalizing the source flux to 10~ 8 </>_s ph(> 100 MeV) cm -2 s" 1 , we have 
<f> a (E,t) = ku~ a v h , where k = lO- 8 a0_ 8 /E loo . For E > 100 MeV, 

S = — / 7 (A o . 2 t 4 )a0_ 8 n-- a84 , (3) 



2 



where g = g{a v ) = 1 — (a^/0.84). Here we use the following conventions 
to define the photon number index a p h, energy index a, and vF v index a u : 
oc — a ph — 1 an d a u = 1 — a = 2 — a p h- For an intrinsic source spectrum with 
photon index a p h = 2, most of the photons are detected near the lower energy- 
threshold. 



3 Background Counts 

The number of background photons with energies > E observed during time 
At that lie within solid element AQ that is centered in the direction Q of a 
specified source is 

At oo 

B = J dt J dE' AQ(E') A[E', 6{t), <f>(t)) <S> B {E', Q) , (4) 

E 



where <&b{E',Q) is the diffuse 7-ray background in the source direction. A 
minimum level of background is provided by the diffuse extragalactic 7-ray 
background radiation, given by 

$ X (E) = k x u~ aB , 0.4 < u < 700 (5) 

0), where k x = 1.73 ± 0.08 x 10~ 7 ph (cm^s-sr-MeV)" 1 and a B = 2.10 ± 
0.03. Other sources of background, in particular, the galactic diffuse 7-ray 
background, can also contribute; the > 100 MeV photon flux is ~ twice as 
bright at |6| = 35° and over 30 times brighter at the Galactic center than 
at the poles (0). A better estimate of blazar flaring rates must consider the 
galactic diffuse emission. 

Considering only the extragalactic diffuse radiation and letting AQ = AOq, 
the number of background counts detected by GLAST is 

B = IAX0.2UU- 2 27 . (6) 



This expression is valid for 1 ^ u £ 100. 
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4 Five Count and 5a Criteria 



The detection of at least 5 counts implies from equation (4) that S > 5S 5 with 
S 5 > 1 (the 5 count criterion), so that 

S4 

(X o . 2 t 4 )0-8 £ gS 5 u°- S4 - a » . (7) 



This expression fails for hard spectrum sources with ot u £ 0.84, where cutoffs 
or breaks in the higher energy range of the spectrum and high-energy correc- 
tions to the GLAST effective area determine the number of detected source 
counts. Otherwise, most of the source counts are detected at the lower energy 
range near E = E 100 . At energies below 100 MeV, the effective area declines 
more rapidly, so that most photons will be detected at even lower energies 
than 100 MeV unless a v > 0.63. 

Detection of a source with significance > 5a requires that S/ \fB > 5<r 5 . Using 
expressions (4) and (7) for S and B, respectively, gives 

(X o . 2 t 4 )0„ 8 > ^ (^) 2 u-^-V . (8) 



By comparing equations (8) and (9), one sees that detection of soft-spectrum 
(a u £ 0) sources with the LAT at the level of 0„ 8 — 15, corresponding to the 
minimum flux sources that EGRET could detect in a 2 week viewing period, 
are background-limited near 100 MeV (that is, the time-on-source Xq^a is 
larger for the 5a criterion when u ~ 1 than for the 5 count criterion). 

For spectra that are not very soft (i.e., a v £ —0.3), the significance of detection 
increases with energy even as the source counting rate declines. Therefore 
there is a unique energy E = Eiqqu at which both the 5 count and 5a criteria 
for a source emitting at a given flux level 0_§ are satisfied, and this energy 
defines the minimum observing time to satisfy these criteria. For sources with 
ol v £ — 0.3, this energy is 

U = (3^)1/(1.43+0.,) ^ 16 (Q) 

a/ 7 0_ 8 5 5 (0„ 8 5 5 ) - 7 



where the last expression applies in the special case a v = 0, / 7 = 0.68. When 
</>_8 = 15, the energy to meet the two criteria is u = 4, or the greatest 
sensitivity is at 200-300 MeV. For sources softer (more negative) than a u = 
—0.3, energies near 100 MeV give the best sensitivity. 
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Fig. 1. The size distribution of blazar flares with 2- week-average fluxes (</>) measured 
during the phase 1 EGRET all-sky survey ((51). The straight line has slope —3/2. 

The minimum observing time is given by 

(x . 2 u) mm = ^ 9 ^u °- 84 -- - 7 s 5 o - 4 v 5 L18 (i^r L59 , (10) 



where the right-hand-side is evaluated in the special case described above. 
Thus, GLAST will detect high-latitude sources at the level EGRET could 
reach in two weeks in only ~7xl0 4 s = 0.8 days, and it accomplishes this 
task in the scanning mode. Considering a background twice as large, which 
is more realistic for the majority of blazars, gives a detection of 0_s — 15 
sources in « 1.2 days. Consequently, GLAST will scan the entire sky to the 
level of 0_g = 15 in a time slight longer than one day. For this, we can take 
the EGRET phase 1 catalog as a guide to what GLAST will observe. 



5 Comparison with EGRET Phase 1 Catalog 



Phase 1 of the Compton Gamma Ray Observatory observing program took 
place during a 596 day interval from 1991 April 22 - 1992 November 17 (0). 
Fig. 1 shows the size distribution of the two-week average fluxes of radio-loud 
quasars and BL Lac objects measured with EGRET during Phase 1, which 
consists of 44 high-confidence detections from 25 different blazars. Because 
the CGRO Phase 1 was an all-sky survey, it will be less biased to reveal the 
outcome of GLAST scanning observations than later EGRET catalogs where 
the number of detections is biased by the many pointings in the direction to 
previously known bright gamma-ray sources. Even so, the Phase 1 exposure 



5 



10° 
10 4 

10 3 

* io 2 

o 

"in 10 1 



X 



10 
10" 

10 
10 



N 

N 

< N: s 



1000 



-= 100 



- 10 u 



-2 
-3 

0.1 100 1000 



Fig. 2. Minimum required observing time for GLAST in scanning mode to satisfy 5 
count detection, 5 sigma sensitivity criteria for a high-latitude source at the integral 
photon flux level 10~V-8 ph(> 100 MeV) cm~ 2 s _1 . Right axis is u = E/E 100 , the 
photon energy in units of 100 MeV that satisfies these criteria in the shortest time. 
These results are for sources with flat vF v spectra, i.e., a v = 0, and assumes that 
only the extragalactic diffuse background contributes to background counts. The 
solution for u is valid for 1 ;$ u ?S 100. 

map of EGRET is far from uniform. One must keep in mind that the two-week 
fluxes over shorter, brighter flaring phases, so that the estimated blazar flaring 
rate could be increased if the duty cycle of flares was better known citewalOO. 

The EGRET field-of-view was about l/24th of the full sky, so that the Phase 
1 duration was about twice as long as the time required to survey the entire 
sky with two-week pointings. During this time, EGRET detected some 3-10 
sources at the level </>_s ~ 100. Thus we expect GLAST to detect 2-5 such 
bright sources during a one-day scan. Although uncertain in view of the small 
numbers of blazars detected at bright flux levels with EGRET, one or two 
blazars brighter than 0_ 8 £ 200 should be detected with GLAST every 3-4 
days. A comparable number of unidentified 7-ray sources should also reach 
these levels (M). 



6 Discussion 



Fig. 2 shows the amount of time in the scanning mode required for GLAST 
to detect sources at given flux levels, using equations (10) and (11). Once 
the source flux 0_ 8 £ 200, dozens of photon counts can be expected on a 10 4 
second timescale. Changing to pointing mode for a non-Earth-occulted source 
will increase the source counts by a factor of ~ 3-5, with the larger number 
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applying only to those few sources that are favorably located so that Earth 
occultation is not severe. At these flux levels, there are now sufficiently many 
photons to perform time-resolved spectroscopy. Sources exceeding the level 
of 0_§ ~ 200 must be carefully and quickly reviewed to determine whether 
to issue an alert. Coordinated sensitive radio, optical, UV, and X-ray mea- 
surements would greatly strengthen tests of dominant radiation processes and 
inferred locations of the 7-ray emission sites. At these bright flux levels, 77 
absorption and spectral evolution can be studied to deduce Doppler factors 
and monitor Doppler factor variations 

Fig. 2 also shows that GLAST in the scanning mode will be sensitive to high- 
latitude 7-ray sources at the level of 0_§ = 2 and 0.2 over w 2 weeks and 2 years 
observing time, respectively. Thus GLAST will reach flux levels of ~ 6 x 10 -13 
ergs cm -2 s _1 , fifty times fainter than EGRET's limiting sensitivity. As can 
be seen from Fig. 2, good low-energy (< 100 MeV) response and calibration 
of the LAT on GLAST are important for identifying transients on the shortest 
time scales. 
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